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ABSTRACT 



Ultra Steep Spectrum (USS) radio sources have been successfully used to se- 
lect powerful radio sources at high redshifts (z > 2). Typically restricted to 
large-sky surveys and relatively bright radio flux densities, it has gradually be- 
come possible to extend the USS search to sub-mJy levels, thanks to the recent 
appearance of sensitive low-frequency radio facilities. Here a first detailed analy- 
sis of the nature of the faintest USS sources is presented. By using Giant Metre- 
wave Radio Telescope and Very Large Array radio observations of the Lockman 
Hole at 610 MHz and 1.4 GHz, a sample of 58 USS sources, with 610 MHz inte- 
grated fluxes above 100 fiJy, is assembled. Deep infrared data at 3.6 and 4.5 /im 
from the Spitzer Extragalactic Representative Volume Survey (SERVS) is used 
to reliably identify counterparts for 48 (83%) of these sources, showing an aver- 
age total magnitude of [3.6]ab = 19.8 mag. Spectroscopic redshifts for 14 USS 
sources, together with photometric redshift estimates, improved by the use of the 
deep SERVS data, for a further 19 objects, show redshifts ranging from z = 0.1 
to z = 2.8, peaking at z ~ 0.6 and tailing off at high redshifts. The remaining 
25 USS sources, with no redshift estimate, include the faintest [3.6] magnitudes, 
with 10 sources undetected at 3.6 and 4.5 /xm (typically [3.6] > 22 — 23 mag, 
from local measurements), which suggests the likely existence of higher redshifts 
among the sub-mJy USS population. The comparison with the Square Kilometre 
Array Design Studies Simulated Skies models indicate that Fanaroff- Riley type I 
radio sources and radio-quiet Active Galactic Nuclei may constitute the bulk of 
the faintest USS population, and raises the possibility that the high efficiency of 
the USS technique for the selection of high redshift sources remains even at the 
sub-mJy level. 

Subject headings: surveys — galaxies: active — galaxies: evolution — galaxies: 
high-redshift — radio continuum: galaxies 



INTRODUCTION 



High-redshift radio galaxies (HzRGs) are amongst the most luminous galaxies and seem 
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). The presence of such huge amounts of star 



formation at high redshift, and the fact that we know that powerful radio galaxies in the 



lower redshift Universe ar e hosted by massive ellipticals (IBest et al.lll998t iMcLure fc Dunlop 



Sadler 


2003; 


Dunlop 


2003) 



used as beacons to ident ify over-densities in the distant Universe, i .e. proto-cluster environ 



ments at z ~ 2 — 5 (e.g. IStevens et al.ll2003l ; IVenemans et al.l 120071 ). Identifying and tracing 



the evolution of HzRGs thus offers a unique path to study galaxy and large scale structure 
formation and evolution from the earliest epochs. 

One of the most successful tracers of HzRGs relies on the relation between the steep- 
ness of the radio spectra and redshift (e.g. iTielens et al.lll979l ; IChambers et al.lll996l ). Al- 
though an ultra-steep (radio) spectrum (USS; a < — 1 with S oc u a ) does not guaran- 
tee a high redshift, and the USS selection may actua lly miss a large fraction of HzRGs 



Waddington et al.l Il999l ; iJarvis et al.l l200ld 120091 ) , a higher fraction of high-redshift 



(e.g. 

sources can be found among those with the steepest radio spectra, a nd most of the radio 



galaxies known at z > 3.5 have been found us in g the USS criterion (IBlundell et al 



De Breuck et al. 1998 



Cruz et al 



2006 



2000 



Milev fc De Breuck 



Jarvis et al. 2001a. b; De Breuck et al. 2002a; Jarvis et al. 



1998; 
2004; 



20081 ). A satisfactory explanation for this trend is still 
missing. One possibility is that the steeper radio spectral index arises from a combination of 
an in creased spectral curva ture with redshift, and the redshifting of a concave radio spectrum 



(e.g. 



Krolik fc Chenlll99lf ) to lower radio frequencies. Another explanation, increasingly be- 
ing supported by the observations, suggests that steeper spectral indexes are linked to radio 
jets expanding in dense environments, a situation o ne would more frequently encounter in 



proto-cluster environm ents in the distant Universe (IKlamer et al.l 12006; Br yant et al.l 12009 



Bornancini et al.ll2010l ). A recent study of submillimetre galaxies (llbar et al.ll2010l ). known 
to be high-redshift (z ~ 2) massive objects linked to the galaxy formation process in proto- 
clusters, seems to support this latter hypothesis: while the average slope of their radio contin- 
uum emission is indistinguishable from that of local star-forming galaxies (o^io'mhz 2 = 

-0.75, 

with standard deviation of 0.29), those sources with an active galactic nuclei (AGN) domi- 



nated mid-IR spectrum typically show steeper radio spectral indices (agio° < 



-1.0). 



Until recently, widefield (over several tens or hundreds of square degree) and, conse- 
quently, only moderately deep radio surveys h ave been used to s elect and stu dy the USS 
radio population; for example, the sources of iDe Breuck et al.l (120001 l2002al M 120041 ) do 
not have 1400-MHz flux densities below ~10-15mJy. This raises the question of whether 
a significant population of apparently fainter USS sources either at even higher redshifts 
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and/or with lower AGN power is being missed by these surveys. At sub-mJy levels the 
issue could actually become more complex, as the population mix of radio sources changes 
dramatically from that observed above a few mJy, with an increasingly large contribution 



from the evolving star-forming galaxy populat i on fe.glJarvis fc Rawlings 1120041 ; lAfonso et al. 



20051 ISimpson et aDhood : ISmolcic et aDbo08l : Iwilman et ap2008h . With the recent appear- 



ance of sensitive low-frequency radio facilities such as the Giant Metrewave Radio Telescope 
(GMRT), and with LOFAR being commissioned, it has becom e viable to extend the search for 



and subsequent study of USS sources to the /xJy regime (e.g. iBondi et al.l 120071 ; lOwen et al. 



2009; 


Afonso et al. 


2009; 


Ibar et al. 


2009) 



but still lacking, as the deepest possible supporting multi-wavelength data are an absolute 
necessity. 

The infrared (IR) regime is, in this respect, of paramount importance: not only are 
USS sources potentially at very high redshift, but they are also potentially dusty, which 
means they will be increasingly difficult to detect at optical wavelengths. Wide-area deep IR 
surveys, which have became possible over the last few years with the Spitzer Space Telescope, 
are now allowing to significantly improve the study of USS sources. Of particular relevance is 
the recent transition of Spitzer to the post-cryogenic, or "warm" mission, phase. This change 
has allowed much larger surveys with the two shortest wavelength channels (3.6 and 4.5 //m) 
of the Infrared Array Camera (IRAC), than was sensible during the cryogeni c portion of the 



missi on. The Spitzer Extragalactic Representative Volume Survey (SERVS - iMauduit et al. 
201 ll ) is such a wide-area survey, covering 18 deg 2 to w 1—2 /xJy in both the [3.6] and 



[4.5] bands, over five of the best-studied extragalactic fields: the European Large Area ISO 
Survey fields Nl and SI, the XMM-Newton Large Scale Structure survey, the Chandra Deep 
Field South and the Lockman Hole. 

All these fields were, or are being covered at radio wavelengths to very deep levels (see 
Mauduit et al.ll201lL for details). We have selected one of them, the Lockman Hole, for a first 
detailed analysis of the faintest USS sources, given the availability of uniform coverage at two 
radio wavelengths, 1.4 GHz and 610 MHz, the early availability of the entire set of SERVS 
data in this field and the existence of a large body of supporting multiwavelength data. This 
paper is organised as follows. In Section [2] we describe the radio data and the USS sample 
selection. The multiwavelength data considered in this work is detailed in Section [31 while in 
Section|l]the identification of USS sources at IR wavelengths is presented together with some 
indications of the nature of these sources from IR colour-colour diagnostic plots. The results 
of spectroscopic redshift determinations and photometric redshift estimates are considered in 
Section 0, along with the implications for USS source characterization. Section [6] compares 
the results from this work with SKADS Simulated Skies radio source population models and 
addresses the efficiency of the USS criteria at sub-mJy radio fluxes for the detection of very 
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high-redshift sources. Finally, our main conclusions are summarized in Section [7J 

Throughout this paper we adopt H = 71 kms -1 Mpc -1 , Q M = 0.27, and = 0.73. 
All quoted magnitudes are in the AB system. 



2. USS sample selection 



The Lockman Hole radio coverage at 610 MHz and 1.4 GHz from llbar et al.l ( 120091 ) was 
chosen for the selection of USS radio sources. Besides the depth and resolution of these 
surveys (reaching an r.m.s. of 6/iJybeam~ 1 at 1.4 GHz, with a restoring beam of 4.3" x 
4.2" and 15/iJybeam -1 at 610 MHz, with a 7.1" x 6.5" restoring beam), they have been 
reduced in a similar way, using the same procedures, and with source catalogues which are 
easily comparable to each other. Problems associated with matching sources observed at 
different resolutions at different frequencies, for example, or using different source detection 
algorithms, are thus minimised. 

The area covered by both surveys is slightly different: the 1.4-GHz coverage with the 
VLA amounts to ~0.56deg 2 while at 610 MHz the GMRT coverage is nearly ldeg 2 . In 
order to select USS sources we require coverage at both frequencie s, and thus restricted 
our search to the smaller 1.4-GHz area. As indicated in llbar et al.l ( 120091 ). the 610-MHz 
catalogue contains 1,236 c omponents with in the area covered by the 1.4-GHz observations. 
After component merging (jlbar et al.ll2009l ) we are left with 1,213 individual radio sources. In 
order to perform a more robust selection of USS sources we further restrict the sample in two 
ways: firstly, to avoid the unreliable spectral indices measurements of the faintest detected 
sources, we limit the 610 MHz radio sample to integrated flux densities of S6iomh z > 100 /xJy; 
secondly, we reject all sources at the edge (within ~ 20 ") of the 1.4-GHz map, where the 
rapidly varying noise can bias the measured 1.4-GHz flux. This selection leaves a total of 
928 (610-MHz selected) radio sources. After cross-correlation with the 1.4-GHz catalogue, 
we find spectral indices, a^io (using S u oc u a , where S u represents the total integrated flux 
density at frequency v) for 662 of these, with an upper limit (non-detection at 1.4 GHz) 
obtained for 193 sources. We discard the remaining 73 sources, as th ey have an unreliable 
radio spectral index determination (flagged entries in llbar et al.l 120091 ) . Figure [1] represents 
the flux density at 610 MHz as a function of spectral index and the spectral index distribution 
for the sources in this field. 

For the definition of the USS source sample we have considered commonly used pa- 
rameters. These are by no means unifor m, not in the frequen cies used nor in the spectral 



index threshold - e.g. aff?° < -0.981 dBlundell et al.l Il998h . ajf < -1.2 (ICohen et al. 
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200 41), eg < -1.3 (iDe Breuck et all 12004]) . a\iT < -1.0 f lCruz et al.ll2006l ). affi < -1.0 



( IBroderick et al.l 120071 ). In this work we have adopted a selection based on a conservative 



value of ag^j < —1.3, which results in a sample of 65 USS sources, 10 of which are not de- 
tected at 1.4 GHz (i.e. only have an upper limit in ctgffj based on a 5<r threshold at 1.4 GHz). 
Given the sensitivity of the USS criteria to the measured radio fluxes, we have visually in- 
spected the radio images to identify and reject cases where the USS classification could 
be questionable. We thus rejected all cases where the comparison of the higher-resolution 
1.4-GHz map with the 610-MHz map indicates that the radio emission at these frequencies 
traces different structures - for example, when the 1.4-GHz map shows a pair of close neigh- 
bours that appear blended in the 610-MHz map. We have found 7 such cases. The final 
USS sample in the Lockman Hole thus consists of 58 sources (including 10 undetected at 
1.4 GHz), details of which are presented in Table [TJ 

Given the flux measurement errors, the relative proximity of the radio frequencies em- 
ployed, and the consequent relatively high spectral index uncertainties, it is important to 
characterize, in a statistical sense, the reliability of the USS sample. In particular, it is 
important to evaluate what is the expected ratio of non-USS sources that scatter into the 
USS regime due to measurement errors. We begin by considering an intrinsic spectral in- 
dex distribution for 610 MHz-selected radio sources described by a normal distribution of 
( a 6io°) i^a = —0.80 ± 0.33 and a distribution for the spectral index measurement errors 
(Aa) described by a normal distribution of (Act) ± a^ a = 0.30 ± 0.12. The latter is drawn 
directly from the observed measurement errors for 610 MHz sub-mJy sources in the Lock- 
man Hole and the former is chosen as to be compatible w ith the observed spectral index 



distribution for the sub-mJy radio sample (see Table 7 in llbar et al.ll2009l ). Under these 
assumptions, we find that a significant number of sources with measured spectral indices 
a lto° — — 1-3 are likely to have intrinsic spectral indices which are higher. This is not a 
surprise, but the expected behaviour from the selection of s ources in the tail of the spectral 



index distribution (see, for example, iDe Breuck et al.ll2004j ). and can only be minimised by 



establishing a conservative cut in a and/or by trying to minimise the error affecting the 
a measurement - by working at strong flux densities and/or by choosing widely separated 
radio frequencies for the measurement of a. For the current work, the availability of quality 
radio data in the Lockman Hole and the objective of studying the faintest USS sources led to 
the choice of a conservative radio spectral index threshold (a^io < —1.3) for their selection. 
From the simulations described above, we expect only about a third of the sample in Ta- 
ble [T]to have an intrinsic radio spectral index below or equal to —1.3. However, most of the 
other sources can still be considered USS sources: around 75% of the sample is expected to 
have an intrinsic radio spectral index of a^ < —1.0, still a commonly used USS selection 
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thresholcQ. This justifies our earlier conservative choice for the a^g threshold, and indicates 
that the USS sample in Tabled] is appropriate for the statistical study of the faint radio USS 
source population - keeping in mind that around 14 sources could have celf® > —1.0. In 
order to help identifying any bias from the expected fraction of non-USS sources (those with 
an intrinsic radio spectral index > —1.0) included in the sample, we define a control 

sub-sample - labeled as "A" in Tableland henceforward designated USS(A) subsample - 
including only the 30 sources with the most secure USS classification. 

We note that the assembled USS sample is not intended to be complete: the uncertain- 
ties in the spectral index determination will also scatter real USS sources into a non-USS 
classification. Completeness is nevertheless a requirement beyond the scope of this work. 



3. Optical and Infrared multi- wavelength data 



The Lockman Hole is one of the most extensively observed regions in the sky. Several 
surveys have obtained deep data on this field, from the X-rays to the radio. In this work, 
we make use of various multi-wavelength datasets for the identification of USS faint radio 
sources and the estimate of their redshift distribution. For the identification we use the IR 



data from SERVS ( iMauduit et al.l l201lh . which provides deep and uniform coverage of the 
entire region at 3.6 and 4.5 fim. For the redshifts of these sources we use the spectroscopic 
data available on this field and complement it with photometric redshift estimates for sources 
with five or more broad-band photometric detections in the optical to IR (0.36 — 4.5 /zm) 
wavelength range. For this we use the ground-based optical imaging obtained with the Wide 
Field Camera at the Isaac Newton Telescope (Uq 'r'i'Z bands), to 5a point sou rce depths of 
24.1, 24.9, 24.3, 23.6, 22.1 mag (AB) respectively fiGonzalez-Solares et al.lhoioh . the ground- 
based near-IR data, a t J and K bands, from the UKIDSS Deep Extragalactic Survey (DXS - 
Lawrence et al. I l2007h reaching a limiting magnitude of Kab ~ 21.3 mag and the IR data from 
SERVS. Finally, we also use the Spitzer Wide-area Infrared Extragalactic survey (SWIRE - 
Lonsdale et al.ll2003f ) measurements at wavelengths of 5.8 and 8.0 /zm, to explore the infrared 
colours of the USS sources. This vast body of wide-area multi-wavelength data was merged 
into a S ERVS-selected mult i- wavelength catalog, or SERVS Data Fusion, using a 1.5" search 
radius (IVaccari et al.l 120101 ) . Deeper optical and near-IR observations of the Lockman Hole 



do ex ist, for parts of the field, e.g. Subaru i maging in the r, i and z bands ( jSzokoly et al. 



20101 ). and LBT observations at U, B and V (IRovilos et al.ll2009l ). These will be used for the 



1 We note that even increasing the spectral index measurement errors to (Aa) ± a& a — 0.40 ± 0.12 would 
still result in about two thirds of the selected sources having intrinsic radio spectral index of a|fo° < —1-0. 
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study of specific USS sources in future papers, but the statistical study of the USS sample, 
collected over a wide-field, requires a more homogeneous approach which is facilitated by 
the wide-area datasets adopted here. 



4. The IR counterparts of USS faint radio sources 
4.1. Identification rate 

Counterparts for the samp le of USS sources were searched for in the SERVS images using 



the likelihood ratio method of ISutherland fe Saunders! (119921 ) . Due to the better resolution 



delivered at the higher frequency, the nominal radio position was taken as the 1.4-GHz 
position, whenever a detection at this frequency was available (with the 610- MHz position 



adopt ed otherwise). The radio positional uncertainty was estimated following llvison et al. 



( 120071 ). We considered further that 3cr position al uncertainties of less than 1" are unrealistic 



given the possibility of systematic offsets (e.g., iMorrison et al.ll2010l). For each USS source. 



the 3 .6 Atm band identification with the highest reliability (1Z, see ISutherland fe Saunders 



19921 ) above 75% was taken as the real counterpart. The 4.5 /im band image was inspected 
for identifications not present in the shorter wavelength band, but none were found. All 
identifications were inspected visually to check for special situations in which the likelihood 
ratio method would not apply, as in the case of non-independent sources (either in the radio 
or in the IR). 

Using the deep SERVS observations 48 out of the 58 (83%) USS sources possess an 
IRAC identification with 1Z > 75%. The median AB magnitude for the USS source sample is 
[3.6] = 19.7mag (an average value of [3.6] = 19.8 mag is found among detected counterparts). 
For the ten USS sources with no IRAC identification an estimate of the 5 a [3.6]-magnitude 
lower limit has been estimated by using several measurements around the radio position. 
These results do not seem to be influenced by the contamination of the sample by non-USS 
sources. If we consider the more robust USS (A) sub-sample, we find an equally high IRAC 
identification rate of 80%, and a very similar median magnitude for the counterparts of [3.6] = 
19.6 mag (also an average of [3.6] = 19.6 mag considering only the detected counterparts). 

The identification rate for the USS sample using the infrared data in SERVS is compa- 
rable to that for the overall radio population in the Lockman Hole: 773 out of the 848 (91 per 
cent) sources with a radio spectral index determination or upper limit have an IRAC identi- 
fication. It is important to note that chance coincidences have a limited impact in these high 
identification rates, as we have rejected all counterpart candidates with reliabilities below 
75 per cent (and the vast majority of the identifications show very high reliabilities, 1Z > 95 
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per cent). The false identification rate should thus be at the few per cent level. We have 
confirmed this result by running identification searches for random radio positions, obtained 
from the real USS sources positions (and assuming the same error regions) by adding a 
random distance between 30 and 45" along a random direction. The simulations returned 
random identification rates of between 2 and 4 per cent, as expected. 

The identification rate is also surprising in another way: when different radio spec- 
tral index ranges are considered, the identification rate always remains at the ~90-per-cent 
level (Table [2]). This result is in contrast to previous work with USS sources at higher ra- 
dio flux densities and brighter optical/IR magnitu de limits, showing an identification rate 



which decreases with steeper radio spectral index ( IWieringa fc Katgertlll99ll ; llntema et al. 



20101 ). Such indications are usually taken to imply that sources with steeper radio spec- 
tral indices are, on average, at higher redshifts. Our high identification rate for the sub- 
mJy ra dio USS population does not, however, imply otherwise. The situation is akin to 



that of iDe Breuck et al.l (J2002aJ) which, by using relatively deep -fT-band observations (to 
K ~ 22 mag), achieve a near-IR identification rate of 94 per cent for their radio bright 
(> lOmJy at 1.4 GHz) USS sources. The dependence of the identification rate vs. radio 
spectral index relation on the optical/IR magnitude limit is clear from Figure [2J When 
the limiting magnitude for the infrared counterparts is set at [3.6] = 21, 20 and 19 mag 
not only does the identification rate show a significant decrease for all spectral indices but, 
by [3.6] = 19 mag, a more pronounced decrease of the identification rate of steeper spec- 
trum sources seems to be settling in (23 per cent for ct\^ < —1.3 but ~ 40 per cent for 
a 6io° > — 0.9). This is a direct result of fainter [3.6] magnitudes for the sub-mJy USS 
population when compared to the flatter (a^ > —0.9) radio spectral index population: 
as mentioned above, the median magnitude for USS (measured c4io° — — 1-3) sources is 
[3.6] = 19.7mag (average of [3.6] = 19.8 mag for detected counterparts) while amongst 
sources with a measured radio spectral index of c^q > —0.9 the median magnitude is 
[3.6] = 19.3 mag (average of [3.6] = 19.2 mag for detected counterparts). In parallel to the 
"traditional" (above the mJy level) USS population, this is also compatible with higher 
average redshifts among the steepest spectrum sources. 

The issue of the contamination of the USS sample by non-USS sources should also be 
considered here, as the spectral index measurement errors could affect the observed behaviour 
of the identification rate with spectral index. In general, measurement errors in the spec- 
tral index determination would be expected to flatten out any intrinsic correlation between 
identification rate and spectral index, if one exists, given the more significant scattering of 
presumably higher-identification rate non-USS sources (which exist in higher relative num- 
bers) into the low-identification rate (and less populated) USS regime. The question then 
becomes how much of the almost flat behaviour of the identification rate with spectral index 
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present in Figure |2] is due to contamination of the USS sample by non USS sources. This 
can be tackled by simulating the effects of applying the distribution of spectral index errors 
in the current radio sample (considering only the sub-mJy regime, being the relevant one) to 
some intrinsic identification rate vs. radio spectral index relation. In Figure |3] we show the 
intrinsic identification rate vs. radio spectral index relation which, for an intrinsic spectral 
index distribution of 610 MHz-selected radio sources described by a normal distribution of 
( a lw°) =t °a — —0.80 ± 0.33 and a distribution for the spectral index measurement errors 
(Act) given by a normal distribution of (Aa) ± a^ a = 0.30 ± 0.12 (see section [2]), would 
result in the observed relation. Steeper identification rate vs radio spectral index relations 
have been explored but found incompatible with the observations. The effect of the spectral 
index measurement errors or, equivalently, the effect of the contamination of the USS regime 
by non-USS sources is, under these assumptions and for the purpose of the identification 
rate vs radio spectral index relation, virtually inexistent, and cannot be the reason for the 
flat curve in Figure EJ As detailed above, this should then be mostly a result of the depth 
of the infrared observations employed here. 



4.2. IR colour diagnostics for USS sources 



Given the Spitzer coverage of the Lockman Hole region, both with SERVS and SWIRE, 
IR colours are available for many of the USS sources. This allows for the use of IR colour 
diagnostics plots, which may offer further clues about the nature of these sources. Figure 0] 
shows the location of radio sources in the Lockman Hol e with appropriate IR detections in 
two of the most commonly used IR diagnostic diag rams (ILacv et al.ll2004i IStern et al.ll2005l ) 
together with the recent KI diagnostic diagram of iMessias et al.l (120 111) . USS sources show 
a wide range of IR colours, with a significant fraction being associated with AGN colours - 
among those USS sources with simultaneous detections at the relevant infrared bands about 
two-thirds fall in the AGN region for the Lacy et al. criteria, while around a third do so 
for the Stern et al. and the Messias et al. criteria. In spite of the small numbers of sources 
with appropriate photometry, the USS (A) subsample does not appear to differ significantly 
from the full USS sample. While these criteria are known to be contaminated by non-AGN 
galaxies, in particular at high (z > 2.5) redshifts, a nd are certainly not complete (i.e., AGN 
sources do appear outside the designated regions, see Donley et all2008~ ; Barmby et al. 2008 ). 
they provide a strong indication that the USS population contains a substantial amount of 
IR-detected AGN. We note, nevertheless, that different accretion modes may be found in 
radio selected AGN, and in partic ular a radiatiyely inefficient "radio mode" may not be 



detectable at IR wavelength s (e.g., ICroton et al.l l2006t lHardcastle et al.ll2007t iTasse et al. 
20081 : iGriffith fc Stern! l^oL 
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Redshift distribution 



5.1. Spectroscopic redshifts 



Several spectroscopic o bservations have targ e ted specific classes of sources in the Lock- 



man Hole in the past (e.g. iLehmann et al.ll200ll ; IChapman et al.ll2005l . for X-ray and sub- 



millimetre sources, respectively), but, besides the SDSS coverage of the brighter optical 
sources, only recently extensive deep spectroscopic surveys have been initiated. Among the 
58 USS sources in this work, we have found two spectroscopic redshift deter minations from 



publi s hed follow-up observa tions of X-ray sources in the ROSAT Deep Survey (ISchmidt et al. 



19981 ; ILehmann et al.ll200l|). six from an ongoing MMT-hectospec spectroscopic survey of 
SWIRE 24 [im sources (IHuang et al.ll2010l ). and four were found among a recent Keck spec- 
troscopic survey of faint radio sources in this field (P.I. Scott Chapman). Finally, two sources 
had recent mid-infrared spectroscopy observations obtained with the Spitzer-lnfvaRed Spec- 
trograph (IRS ) , as p art of an analysis of candidate AGN-dominated sub-millimetre galaxies 
( ICoppin et al.ll2010f ). These observations originate the two current highest spectroscopic 
redshifts in the USS sample, at z — 2.56 and z = 2.76. Table |3] shows details for these 14 



sources. 



5.2. Photometric redshifts 



For the remainder of the USS sample it is possible to obtain photometric redshift es- 
timates. So me of the sources have had a photometric redshift estimate from the SWIRE 
survey, from iRowan- Robinson et al.l (120081 ). The procedure there employs both galaxy and 
quasi-stellar object templates applied to data at 0.36 — 4.5 /im and a set of four IR emission 
templates fitted to IR-excess data beyond 3.6 fim. We find a total of 12 objects with a pho- 
tometric redshift derived using 5 or more photometric bands, with 5 objects overlapping the 
current spectroscopic redshift sample. H owever, near-IR data, u n availa ble at the time for 
the Lockman Hole, could not be used by IRowan- Robinson et al.l (120081 ) . Given the recent 
availability of data from the UKIDSS Deep Extragalactic Survey on the Lockman Hole and 
the ne w, deeper, SERVS dat a, we have inc luded these photomet ric bands and used both 
EAZY ( Brammer et al. 2008) an d HyperZ (jBolzonella et al. 2000 , following the procedure 



Maraston et al 



20061 ) to find photometric redshifts estimates for a larger number 



outlined in_ 

of sources. We use aperture corrected magnitudes measured in a radius of 1.6", 2.0" and 
1.9" for INT WFC, UKIDSS and IRAC respectively. The strength of EAZY's template error 
function was adjusted to provide the best fit to the existing spectroscopic data; although few 
objects have reliable spectroscopic redshifts, the chosen value was very similar to that used 
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by lSimpson et al.l ( 120101 ) for the 100-/iJy 1.4 GHz radio sample of ISimpson et al.l ( 120061 ) 



Accepting only those photometric redshifts resulting from 5 or more photometric mea- 
surements in the 0.36 — 4.5 /zm range, we obtain estimates for 30 sources, 11 of which have 
an available spectroscopic redshift measurement. The comparison between photometric and 
spectroscopic redshifts is indicated in Figure |5j and shows a good agreement out to z ~ 2. 
The performances of EAZY and Hyper Z are similar, w ith a nma d (using the normalized median 



absolute deviation as defined in lBrammer et al.ll2008l ) of |Az|/(l + 2: S p ec ) of 0.06 and 0.07, re- 



spectively. There is one "catastrophic" outlier (|Az|/(l + z spec ) > 5cr nmad ) for the photomet- 
ric redshift of LH610MHzJ105407.0+573308 as given by EAZY (|Az|/(l+z spec ) = 6.3 o nmad ). 
This source has the second highest spectroscopic redshift of the USS sample, z spec = 2.56, a 
value obtained from IRS mid-infrared spectroscopy. A detailed inspection of the result- 
ing photometric redshifts probability distributions for this source shows broad and flat 
distributions between z ~ 1 and z ~ 3, with high probability secondary solutions at 
z ~ 2.3 — 2.4. This highlights the unavoidable limitations of photometric redshift esti- 
mates for some sources, in particular at higher redshifts. In Table [1] and in what follows we 
adopt the HyperZ estimates, noting there would be no significant change by the adoption of 
the photometric redshifts from EAZY. 

The redshift distribution for the sub-sample of USS sources with a redshift determination 
is presented in Figure [61 It shows a significant presence of high redshift sources, with a peak 
at z ~ 0.6 and extending beyond z ~ 2. The redshift distribution appears to be unchanged 
if one considers only the USS(A) subsample (inset in Figure ED, which points to a limited 
effect of the small number of real non-USS sources expected to be included in this sample. 



The 



J 610 MHz 



radio monochromatic 
. lO^WHz" 1 and L 



610 MH 



luminosities for these sources vary between 
- 10 26 WHz -1 : while the lower values could be 



explained by lower power AGN or even star-forming galaxies (e.g. lAfonso et al.ll2005l their 
Figure 7) - although the latter would pose problems for the interpretation of the steepness 
of the radio spectrum - the higher values are clear indications of powerful AGN activity, 
already above the FRI/FRII luminosity break. Although beyond the scope of the present 
paper we should stress other clear indications for the existence of powerful AGN in the USS 
sample - two of the spectroscopic red shifts were obtained by follow-up observatio ns of X-ray 
sources in the ROSAT Deep Survey (ISchmidt et al.lll998l ; iLehmann et al.ll200ll ). and show 
broad optical lines which classify them as type I AGN; also, as mentioned above, two other 
sources have been observed by Spitzer-IKS as candidate AGN-dominated sub-millimetre 
galaxies, and SED fitting to their mid-in frared emission inde ed indicates substantial AGN 
contributions over this wavelength range (jCoppin et al.ll2010l ). 



One should recall that a redshift estimate exists for only 33 of the 58 USS sources, 
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leaving more than 40 per cent of the sample with no redshift information. Figure [7] shows 
that these objects correspond to the faintest [3.6] magnitudes, being potentially the highest 
redshift sources in our sample^. In particular, the 10 USS sources with no IRAC detection are 
particularly interesting candidates for very high redshift objects. For these sources, the ratio 
of radio to 3.6 jum fluxes can easily be above 50. Such USS sources are then also classified as 
Infrared- Faint Radio Sources ( IFRS), a rare clas s of object first identified in the Australia 
Telescope Large Area Survey (INorris et al.l 120061 ) . Recent follow-up studies have provided 
strong indications that the IF RS population probably contains a majority of obscured, high 



redshift, radio- lou d galaxies (INorris et al. 



2007; 



Middelberg et all 120081: iHuvnh et al 



2010 



Norris et al.ll2010i ). Interestingly, iMiddelberg et al.l (120101 ) . by studying the radio spectral 



indices of IFRS, also finds a significant overlap with the USS population (6 out of 17 IFRS 
were found to be USS sources). This argues for the likelihood of high redshifts for the 10 
USS sources in the Lockman Hole without IRAC-SERVS counterparts. 



6. Comparison with SKADS Simulated Skies radio population models 

While investigating the nature of USS radio sources at sub-mJy radio flux levels, and 
comparing it with the better known population revealed at brighter radio fluxes (above tens 
or hundreds of mJy), it is of the utmost interest to consider the predictions from radio 
population models. Unfortunately, radio spectral index information is usually not robustly 
included in such models, a result of a severe lack in understanding the physical origin of radio 
spectral indexes that differ from the nominal a ~ —0.7 synchrotron value. Keeping such 
limitation in mind, we use the Square Kilom etre Array Design Study (SKADS) Simulated 



Skies (S 3 ) simulations (IWilman et al.l 120081 ) for two distinct situations: (a) Survey Model 
I - a radio survey reaching a detection sensitivity of 100 fiJy at 610 MHz over 0.6 square 
degree, which will be directly comparable with the current work; (b) Survey Model II - a 
radio survey reaching a detection sensitivity of lOmJy at 610 MHz over 50 square degree - 
a situation more akin to the typical USS source search above the tens of mJy level. The 
simulations predict the number of sources, and respective redshifts, revealed by such surveys, 
among the following source types: radio-quiet (RQ) AGN, radio-loud AGN (including both 
FRI and FRII types), Gigahertz Peaked Sources (GPS) and star- forming galaxies. 

In Figure [SJ the redshift distributions for sources detected in both Survey Models are 
represented. The major obvious difference is the source population which the USS criteria 



2 One should also note from Figure [7] that the [3.6]-band distribution of USS sources does not show any 
significant difference to that of the more restricted USS(A) subsample. 



-14- 



will act upon. At high flux densities, the major contributors to the source counts will be 
powerful AGN of the FRI and FRII types, the latter dominating above z ~ 2. When 
sampling sub-mJy radio flux densities, star-forming galaxies along with lower-luminosity 
AGN dominate the source detections, with FRI and RQ sources contributing significantly 
at z < 2. At higher redshifts (z > 2), FRIs dominate the survey model detections, however 
this population is very poorly constrained at such redshifts. 

The S 3 models are compatible with the current work, and also with previous investi- 
gations at higher radio flux densities. If indeed the USS criteria is sensitive to radio jets in 
high density environments, it will pick up a fraction of the powerful AGN population. The 
environment conditions, along with the larger number of powerful FRII sources at higher 
redshifts, would result in a relatively high efficiency of selection of HzRGs in surveys sam- 
pling high radio flux densities (above tens of mJy). On the other hand, at sub-mJy radio 
flux densities, the USS criteria will essentially only be able to select sources from the FRI 
and radio-quiet AGN population, as FRIIs are extremely rare (only ~ 3 — 4 are predicted 
to be detected in the current work, irrespective of radio spectral index). The selection of 
pure star-forming galaxies, increasingly abundant at these faint radio flux densities, seems 
unlikely, given the synchrotron radio spectrum (a ~ —0.7) they usually display. In fact, the 
redshift distribution of the USS sources in the current work (see Figure [6J also reproduced 
in the appropriate panel of Figure [8]) is remarkably similar to the FRI and RQ redshift dis- 
tribution from the S 3 simulations in Figure [8] (Survey Model I), peaking just below z ~ 1 
and tailing-off above this redshift. The ~ 40 per cent USS sources in the Lockman Hole 
sample with no redshift estimates, if indeed at higher redshifts, would match well the model 
predictions, which place several tens of FRI sources at z > 2 for this survey. Interestingly, 
the efficiency of the USS technique at sub-mJy flux densities for the selection of very high 
redshift galaxies may still be significantly high - possibly even higher than that for radio 
surveys at much higher flux densities. Both Survey Models indicate similar numbers of 
detectable powerful AGN at z > 2: around 80 FRIIs for Survey Model II and the same 
number of FRIs for Survey Model I. This raises the expectations for the application of the 
USS technique in the upcoming deepest LOFAR surveys, as it indicates the possibility of 
complementing the sampling of the high-redshift highest power radio sources (FRIIs) at high 
radio flux densities with the selection of co-eval lower power AGN population (FRIs) at sub- 
mJy levels. One should stress, however, that the uncertainties in the SKADS simulations, 
in particular for the population of FRIs at high redshifts, prevent us from drawing definite 
conclusions. Nevertheless, the importance of exploring the sub-mJy USS population with 
the faintest infrared magnitudes is obvious. 
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7. Conclusions 

We present a sample of 58 Ultra Steep Spectrum faint (5*610 mh z > 100 /xJy) radio sources 
in the Lockman Hole. High reliability IRAC-SERVS 3.6 /xm counterparts are found for 48 
(83%) of these sources, with an average [3.6] total magnitude of 19.8 mag (AB). Infrared 
colour-colour diagnostics indicate a significant fraction of IR-detected AGN among the radio- 
faint USS population. Using spectroscopic redshifts for 14 of the USS sources, complemented 
with photometric redshift estimates for a further 19 sources, the redshift distribution of USS 
sources is derived. Redshifts range from z — 0.1 to z — 2.8, peaking at z ~ 0.6 and tailing off 
at high redshifts. The remaining 25 sources include the faintest sources at IR wavelengths, 
and are potentially at higher redshifts. These results are essentially unchanged when the 
sample is restricted to the 30 sources (the USS (A) subsample) with the most secure USS 
classification, a procedure followed here to help identifying any bias from the small fraction 
(an estimated ~ 25%) of non-USS sources scattered to a USS classification due to spectral 
index measurement errors. 

The comparison with the SKADS Simulated Skies models shows an indication that 
FRIs and RQ AGNs may constitute the bulk of the USS population at sub-mJy radio flux 
densities, and points to the possibility that the efficiency of this technique for the selection 
of higher redshift sources may be as high as when applied at much higher radio flux density 
levels. 
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Table 1. Ultra Steep Spectrum Sources in the Lockman Hole 
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10 






1 


005 


LH610MHzJ105101 


9+565750 


A 


688 


121 


—2 ri9 +0 ' 37 

z - uy -0.30 


10:51 


01 


80 


+56:57:50.0 


0.42 


















> 23 


27 














LH610MHzJ105105 


0+571922 


A 


176 


55 


- L40 -o:3o 


10:51 


04 


87 


+57:19:23.2 


0.32 


10 


51 


04 


89 


+57 


19 


23 





19 


88 





22 






1 


795 


LH610MHzJ105108 


4+573344 


A 


374 


97 


-l 62+ ' 27 
1 - DZ -0.25 


10:51 


08 


25 


+57:33:44.9 


0.34 


10 


51 


08 


29 


+57 


33 


45 


2 


19 


36 





45 


1 


540 


1 


645 


LH610MHzJ105108 


9+572534 


A 


136 




< -1.46 


10:51 


08 


91 


+57:25:35.0 


1.45 


10 


51 


08 


80 


+57 


25 


35 


5 


21 


45 


1 


01 










LH610MHzJ105115 


3+571820 


A 


213 


55 


-1 63+ 026 

1.UO_q 2 g 


10:51 


15 


19 


+57:18:20.6 


0.25 


10 


51 


15 


20 


+57 


18 


20 


8 


19 


11 





22 









505 


LH610MHzJ105133 


3+571833 




101 




< -1.32 


10:51 


33 


35 


+57:18:33.4 


0.85 


10 


51 


33 


41 


+57 


18 


33 


7 


21 


91 





56 










LH610MHzJ105134 


5+573218 


A 


181 


54 


-1 46+ ' 36 


10:51 


34 


37 


+57:32:17.2 


0.36 


10 


51 


34 


39 


+57 


32 


17 


4 


19 


67 





29 










LH610MHzJ105146 


8+572032 




155 


45 


-1 48+ - 42 


10:51 


46 


58 


+57:20:32.4 


0.66 


10 


51 


46 


55 


+57 


20 


33 





20 


87 





63 





981 






LH610MHzJ105152 


8+571347 




226 


76 


-1 ,0+ ' 34 
J--JU_ 33 


10:51 


52 


73 


+57:13:47.1 


0.63 


















> 21 


79 














LH610MHzJ105155 


2+570409 




195 


54 


1 - oo -0.43 


10:51 


55 


01 


+57:04:09.6 


0.64 


10 


51 


55 


09 


+57 


04 


09 


7 


19 


33 





65 





421 





405 


LH610MHzJ105158 


0+571940 




159 




< -1.30 


10:51 


58 


03 


+57:19:40.4 


1.28 


10 


51 


57 


77 


+57 


19 


39 


9 


20 


32 


2 


15 


1 


368 


2 


080 


LH610MHzJ105200 


3+571052 


A 


135 




< -1.43 


10:52 


00 


32 


+57:10:52.7 


0.85 


















> 19 


80 














LH610MHzJ105201 


4+565947 


A 


2181 


480 


i oo+0. 30 
-l.»2_ 2g 


10:52 


01 


12 


+56:59:47.6 


2.04 


















> 24 


22 














LH610MHzJ105201 


5+571105 


A 


177 




< -1.64 


10:52 


01 


49 


+57:11:06.0 


0.95 


10 


52 


01 


29 


+57 


11 


08 


1 


18 


90 


2 


67 










LH610MHzJ105207 


8+573346 




170 


57 


, oq+O.39 
1 - OU -0.43 


10:52 


07 


66 


+57:33:46.8 


0.36 


10 


52 


07 


67 


+57 


33 


46 


6 


19 


18 





25 





501 





510 


LH610MHzJ105208 


0+571349 




153 


48 


, OQ+0.31 

i-o»_0.32 


10:52 


07 


91 


+57:13:49.6 


0.28 


















> 22 


62 














LH610MHzJ105208 


1+571321 




140 


46 


1 00+0.43 

1. oo_ 43 


10:52 


08 


08 


+57:13:21.8 


0.65 


10 


52 


08 


17 


+57 


13 


22 


7 


21 


56 


1 


16 











Table 1 — Continued 



Source flag 5 6 10 §1400 a 610° a S u Vos <*SER,VS &SERVS [3.6] ^■[i.t]" 2 spcc z pbot 

(juJy) (/iJy) (J2000) (J2000) " (J2000) (J2000) (AB) " 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



LH610MHzJ105208 


3+571342 


A 


178 




< -1.48 


10 


52 


08 


28 


+57:13 


42 


2 





87 
















> 24 


67 














LH610MHzJ105212 


8+570641 


A 


215 


56 


1 „ 1 +o.48 
- L61 -0.40 


10 


52 


12 


63 


+57:06 


41 


3 





72 


10 


52 


12 


71 


+57:06 


41 


7 


19 


09 





77 









695 


LH610MHzJ105214 


0+571841 




131 


44 


1 OO + 0.44 
- L32 -0.42 


10 


52 


14 


04 


+57:18 


42 








66 


10 


52 


14 


02 


+57:18 


42 


3 


20 


20 





39 










LH610MHzJ105220 


6+573930 


A 


296 


54 


-2-03t£|* 


10 


52 


20 


42 


+57:39 


31 


5 





28 


10 


52 


20 


43 


+57:39 


31 


6 


20 


21 





16 






1 


460 


LH610MHzJ105220 


7+573408 




208 


61 


1 A rr+0 41 


10 


52 


20 


76 


+57:34 


08 


6 





42 


10 


52 


20 


73 


+57:34 


08 


7 


18 


63 





29 






1 


100 


LH610MHzJ105225 


6+571337 


A 


173 


49 


- L52 io:34 


10 


52 


25 


53 


+57:13 


38 


5 





31 


10 


52 


25 


57 


+57:13 


38 


4 


19 


45 





37 





468 





660 


LH610MHzJ105229 


1+574615 


A 


357 


89 


-1 .-rr-r-O 54 


10 


52 


29 


04 


+57:46 


16 


6 





77 


10 


52 


29 


03 


+57:46 


16 


8 


17 


97 





23 









405 


LH610MHzJ105229 


4+574009 




158 


45 


-1.49+^ 


10 


52 


29 


17 


+57:40 


09 


4 





64 


10 


52 


29 


16 


+57:40 


09 


6 


22 


75 





25 










LH610MHzJ105230 


2+570024 


A 


486 


110 


1 rrn + O 54 

-l./»_0.39 


10 


52 


30 


16 


+57:00 


24 


3 





80 


10 


52 


30 


06 


+57:00 


24 


4 


19 


17 





82 










LH610MHzJ105245 


5+573745 


A 


291 


84 


.„ + () 24 


10 


52 


45 


39 


+57:37 


45 


9 





21 


10 


52 


45 


39 


+57:37 


45 


9 


19 


34 





05 


1 


677 


1 


495 


LH610MHzJ105256 


1+574148 




257 


86 


— I.o2_ 30 


10 


52 


56 


07 


+57:41 


48 


6 





31 


10 


52 


56 


04 


+57:41 


47 


7 


19 


13 





97 






1 


795 


LH610MHzJ105257 


8+573058 




115 


33 


Ar .4-Q 54 

- L49 -o:6o 


10 


52 


57 


65 


+57:30 


58 


8 





50 


10 


52 


57 


66 


+57:30 


58 


7 


19 


90 





13 


2 


100 






LiriolOMrizJ lUozoo 


1+5/3616 


A 

A 


246 


/ 1 




10 


52 


58 


11 


+57:36 


16 


9 





39 


10 


52 


58 


19 


+57:36 


1 / 


2 


20 


16 





67 






1 


780 


LH610MHzJ105301 


5+573429 


A 


224 


59 


-1 60+ ' 31 


10 


53 


01 


34 


+57:34 


30 


3 





28 


10 


53 


01 


37 


+57:34 


30 


2 


19 


37 





28 









645 


LH610MHzJ105303 


7+571616 


A 


235 




< -1.68 


10 


53 


03 


65 


+57:16 


17 








65 
















> 22 


93 














LH610MHzJ105308 


4+573557 




132 


40 


'-•^ -0.4.7 


10 


53 


08 


43 


+57:35 


58 


4 





50 


10 


53 


08 


34 


+57:35 


58 


3 


20 


67 





77 










LH610MHzJ105310 


5+572320 


A 


131 


29 


1 81 + ' 46 
-l.»l_0.43 


10 


53 


10 


51 


+57:23 


21 


9 





50 


10 


53 


10 


51 


+57:23 


21 


4 


19 


17 





46 









695 


LH610MHzJ105319 


2+572108 


A 


321 


108 


_1 O1+0.17 
L - ,: "--l).l7 


10 


53 


19 


18 


+57:21 


08 


7 





18 


10 


53 


19 


27 


+57:21 


08 


5 


21 


19 





72 


2 


760 






LH610MHzJ105323 


2+571638 


A 


126 




< -1.35 


10 


53 


23 


16 


+57:16 


38 


4 





74 
















> 22 


94 














LH610MHzJ105332 


0+574145 




264 


79 


_1 44+0-40 


10 


53 


32 


16 


+57:41 


44 


2 





46 


10 


53 


32 


12 


+57:41 


44 


4 


18 


49 





33 





113 





200 


LH610MHzJ105345 


3+572328 


A 


158 


44 


, co+0.40 
— l.OCS_ 42 


10 


53 


45 


18 


+57:23 


29 


7 





36 


10 


53 


45 


18 


+57:23 


29 


3 


20 


81 





39 










LH610MHzJ105354 


5+572340 


A 


154 


40 


-1 62+ ' 46 
1 - D ^-0.44 


10 


53 


54 


31 


+57:23 


40 








42 


10 


53 


54 


38 


+57:23 


39 


9 


19 


41 





59 









705 


LH610MHzJ105358 


1+574001 




293 


95 


-1 36+ ' 39 
1 -°°-0.42 


10 


53 


58 


08 


+57:40 


02 


5 





48 


10 


53 


58 


02 


+57:40 


03 


4 


18 


92 





99 






1 


055 


LH610MHzJ105401 


2+572842 




192 


64 


-1 ,-,+0.35 

1-0 -o.36 


10 


54 


01 


00 


+57:28 


43 


7 





36 


10 


54 


01 


01 


+57:28 


43 


8 


19 


44 





14 









605 


LH610MHzJ105407 


0+573308 


A 


285 


77 


-1 56+ ' 34 


10 


54 


06 


87 


+57:33 


08 


7 





37 


10 


54 


06 


87 


+57:33 


09 


2 


19 


80 





50 


2 


560 


1 


710 


LH610MHzJ105429 


1+572939 




130 


43 


-1 o 9 +0.51 
1 " ; ' z -0.52 


10 


54 


29 


08 


+57:29 


41 


5 





42 


10 


54 


29 


11 


+57:29 


41 


2 


19 


57 





40 






1 


010 


LH610MHzJ105500 


9+573344 




253 


75 


-l.46tg;~ 


10 


55 


01 


06 


+57:33 


45 


9 





50 


10 


55 


01 


02 


+57:33 


45 


9 


19 


89 





36 










LH610MHzJ105502 


3+573224 




297 


98 


-1 oo+0.40 
1 - 3 ' : '-0.39 


10 


55 


02 


27 


+57:32 


25 


9 





44 


10 


55 


02 


29 


+57:32 


26 





18 


80 





17 





421 





450 



Note. — The columns display (1) source name; (2) flag for sources included in a more robust sub-sample "A"; (3) and (4) the flux densities at 610 MHz and 1.4 GHz, 
respectively; (5) radio spectral index and associated error; (6) and (7) radio position from the higher resolution 1.4 GHz data, where available, and 610 MHz otherwise; (8) error 
in the radio position; (9) and (10) the IR position of the SERVS-3.6/im counterpart; (11) 3.6 /im counterpart total magnitude; (12) separation between the radio position and 
the IR counterpart; (13) and (14) the redshift measurement from spectroscopy or photometry measurements, respectively. 



to 

CO 



-24 - 



Table 2. SERVS identification rates of faint radio sources in the Lockman Hole 



1400 
"610 




Number of radio sources 


ID rate 


a > 


-0.1 


40 


98% 


-0.5 < a < 


-0.1 


117 


97% 


-0.9 < a < 


-0.5 


273 


92% 


-1.3 < a < 


-0.9 


177 


91% 


a < 


-1.3 


48 


90% 



Note. — This table only considers sources with a spec- 
tral index determination, i.e., with a 1.4 GHz detection. 
Those sources, and in particular USS sources, with only 
a 610 MHz detection will have a lower radio positional 
accuracy, and, consequently, the probability of chance 
coincidences will be higher. Thus the overall identifica- 
tion rate for the general USS population in this sample, 
including ten 1.4 GHz non-detections, is slightly lower. 
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Table 3. Spectroscopic redshifts for USS sources in the Lockman Hole 



Name 




R 


[3.6] 


z 


spec 




Reference for z spec 






(AB) 


(AB) 










LH610MHzJ104949 


.1+571656 


19 


.61 


18.48 





.351 




Huane et al. ( 20101 


LH610MHzJ105004 


.3+572214 


19 


.80 


18.59 





.244 




Huane et al. ( 20101 


LH610MHzJ105108 


.4+573344 


23 


.11 


19.36 


1 


.540 




Schmidt et al. ( 19981 


LH610MHzJ105146 


.8+572032 






20.87 


0. 


,981 


S. 


Chapman, priv. comm. 


LH610MHzJ105155 


.2+570409 


19 


.96 


19.34 





.421 




Huane et al. f2010l 


LH610MHzJ105158 


.0+571940 


23 


.61 


20.32 


1 


.368 


S. 


Chapman, priv. comm. 


LH610MHzJ105207 


.8+573346 


21 


.71 


19.18 





,501 


s. 


Chapman, priv. comm. 


LH610MHzJ105225 


.6+571337 


21 


.37 


19.45 





.468 




Huane et al. (20101 


LH610MHzJ105245 


.5+573745 


21 


.98 


19.34 


1 


.677 




Lehmann et al. f20011 


LH610MHzJ105257 


.8+573058 






19.90 


2. 


,100 


s. 


Chapman, priv. comm. 


LH610MHzJ105319 


.2+572108 






21.19 


2. 


,76 




CoDDin et al. (20101 


LH610MHzJ105332 


.0+574145 


19 


.55 


18.49 





.113 




Huansr et al. ( 20 101 


LH610MHzJ105407 


.0+573308 


23 


.18 


19.80 


2 


.56 




CoDDin et al. ( 20 101 


LH610MHzJ105502 


.3+573224 


20 


.43 


18.80 





,421 




Huane et al. ( 20 101 



Note. - 



The [3.6] and i?-band magnitudes are total magnitudes. 
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Fig. 1. — Flux density at 610 MHz vs. radio spectral index for faint radio sources in the 
Lockman Hole. Sources are colour coded according to their 1.4 GHz radio flux density (colour 
wedge on the right). Open symbols are upper limits in a, representing sources with no 1.4 
GHz detection. The dashed line represents the adopted lower limit for SmoMHz, while the 
dotted line shows the locus of points with Suoomhz — 30 /iJy, corresponding to the llbar et al. 



(120091 ) 5<T peak flux limit at this frequency. Also shown (top) is the distribution of radio 
spectral indexes, with the hashed histogram representing a upper limits. The median value 
for the spectral index distribution (detections only) is -0.69. 
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Fig. 2. — SERVS 3.6 /xm identification rate for faint radio sources in the Lockman Hole as 
a function of radio spectral index, for different [3.6]-band AB magnitude limits. The error 
bars are indicative of a 5% error in the identification rate in each bin. 
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Fig. 3. — The observed SERVS 3.6 /im identification rate for faint radio sources in the 
Lockman Hole as a function of radio spectral index (solid line) and a simulated intrinsic 
identification rate vs radio spectral index relation (dotted line) that would originate it after 
taking into account the measurement errors for the radio spectral index. This shows that the 
scatering of non-USS sources into the USS regime, due to the measurement errors, cannot 
explain the high identification rate observed for all spectral indices. 
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Fig. 4. — Infrared colour-colour diagnostics for radio source s in th e Lockman Hole. Two 
of the most comm only used are shown, that of ILacy et al.l (12004 ) ( top left) and t hat o f 



Stern et al.l (120051 ) (top right), together with the recent KI diagram of iMessias et al.l ( 120111 ). 



The region delimited by the dashed line is the AGN region. Filled points denote the USS 
sample (outlined the USS (A) subsample), overlaid on the general sub-mJy radio population 
in the Lockman Hole. 
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Fig. 5. — Comparison between spectroscopic redshifts and photometric redshifts obtained 
through three different procedures. See Section [5] for details. 
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Fig. 6. — Redshift distribution for radio-faint USS sources in the Lockman Hole. Filled 
histogram denotes sources with a spectroscopic redshift determination, while the open region 
refers to photometric redshift estimates. A further 25 USS sources (43% of the full sample) 
exist but with no redshift estimate, mostly at fainter [3.6] fluxes and likely to be found at 
higher redshifts. The inset shows the redshift distribution for the USS (A) subsample. 
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Fig. 7. — [3.6]-band AB magnitude distribution for USS sources in the Lockman Hole. 
Sources with no available redshift estimate are indicated by the filled histogram. The filled 
circles and attached arrows indicate the USS sources with no IRAC detection, placed at 
the bin location of their 5 a [3.6] magnitude lower limit. The inset shows the magnitude 
distribution for the USS (A) subsample. 
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Survey Model I: 0.6 deg 2 , S 610 >100 AiJy Survey Model II: 20 deg 2 , S 610 >10 mJy 




Fig. 8. — Predictions from the SKADS Simulated Skies models for the redshift distributions 
of radio source populations, irrespective of their radio spectral indices, for two different 
types of radio survey: Survey Model I (left) is a radio survey reaching a detection sensitivity 
of 100 /iJy at 610MHz over 0.6 square degree, similar to the Lockman Hole radio survey 
considered in the current work; Survey model II (right) is representative of a wider (50 square 
degree) and shallower (lOmJy detection limit) survey, comparable to surveys frequently used 
to search for USS sources with radio flux densities above the tens of mJy level. The plot for 
Survey Model I, on the left, also displays the observed redshift distribution for USS sources 
in this work (open histogram), reproducing Figure [6] for a more straightforward comparison. 



